Electronic shot noise in fractal conductors 



00 

o 
o 



oo 



C. W. Groth, 1 J. Tworzydlo, 2 and C. W. J. Beenakker 1 

1 Instituut-Lorentz, Universiteit Leiden, P.O. Box 9506, 2300 RA Leiden, The Netherlands 
^Institute of Theoretical Physics, University of Warsaw, Hoza 69, 00-681 Warsaw, Poland 

(Dated: February 2008) 

By solving a master equation in the Sierpinski lattice and in a planar random-resistor network, 
we determine the scaling with size L of the shot noise power P due to elastic scattering in a fractal 
conductor. We find a power-law scaling P oc L d i~ 2 ~ a , with an exponent depending on the fractal 
dimension df and the anomalous diffusion exponent a. This is the same scaling as the time-averaged 
current /, which implies that the Fano factor F — P/2eI is scale independent. We obtain a value 
F = 1/3 for anomalous diffusion that is the same as for normal diffusion, even if there is no smallest 
length scale below which the normal diffusion equation holds. The fact that F remains fixed at 1/3 as 
one crosses the percolation threshold in a random-resistor network may explain recent measurements 
of a doping-independent Fano factor in a graphene flake. 
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Diffusion in a medium with a fractal dimension is char- 
acterized by an anomalous scaling with time t of the 
root-mean-squared displacement A. The usual scaling 
for integer dimensionality d is A oc i 1 / 2 , independent of 
d. If the dimensionality df is nonintcger, however, an 
anomalous scaling 



A oc tV(2+«) 



(1) 



with a > may appear. This anomaly was discovered in 
the early 1980's 0, H [1 1 @ and has since been studied 
extensively (see Refs. 0,01 for reviews). Intuitively, the 
slowing down of the diffusion can be understood as aris- 
ing from the presence of obstacles at all length scales — 
characteristic of a selfsimilar fractal geometry. 

A celebrated application of the theory of fractal diffu- 
sion is to the scaling of electrical conduction in random- 
resistor networks (reviewed in Refs. [1, 0). According to 
Ohm's law, the conductance G should scale with the lin- 
ear size L of a d-dimensional network as G oc L d ~ 2 . In a 
fractal dimension the scaling is modified to G oc L d f~ 2 ~ a , 
depending both on the fractal dimensionality df and on 
the anomalous diffusion exponent a. At the percolation 
threshold, the known [f| values for d — 2 are df = 91/48 
and a = 0.87, leading to a scaling G oc L~ 097 . This 
almost inverse-linear scaling of the conductance of a 
planar random-resistor network contrasts with the L- 
independent conductance G oc L predicted by Ohm's 
law in two dimensions. 

All of this body of knowledge applies to classical re- 
sistors, with applications to disordered semiconductors 
and granular metals [l(| EH- The quantum Hall ef- 
fect provides one quantum mechanical realization of a 
random-resistor network (l2| . in a rather special way 
because time-reversal symmetry is broken by the mag- 
netic field. Very recently [l3j], Cheianov, Fal'ko, Alt- 
shuler, and Aleiner announced an altogether different 
quantum realization in zero ma gnet ic field. Following ex- 
perimental [14j and theoretical [15| evidence for electron 
and hole puddles in undoped graphene 16], Cheianov 
et al. modeled this system by a degenerate electron gas 



in a random-resistor network. They analyzed both the 
high-temperature classical resistance, as well as the low- 
temperature quantum corrections, using the anomalous 
scaling laws in a fractal geometry. 

These very recent experimental and theoretical devel- 
opments open up new possibilities to study quantum me- 
chanical aspects of fractal diffusion, both with respect to 
the Pauli exclusion principle and with respect to quan- 
tum interference (which are operative in distinct temper- 
ature regimes). To access the effect of the Pauli princi- 
ple one needs to go beyond the time-averaged current I 
(studied by Cheianov et al. 13]), and consider the time- 
dependent fluctuations SI(t) of the current in response 
to a time-independent applied voltage V. These fluc- 
tuations exist because of the granularity of the electron 
charge, hence their name "shot noise" (for reviews, see 
Refs. 0,[]1). Shot 

noise is quantified by the noise power 



P 



dt(6I(0)6I(t)) 



(2) 



and by the Fano factor F = P/2eI. The Pauli principle 
enforces F < 1, meaning that the noise power is smaller 
than the Poisson value 2el — which is the expected value 
for independent particles (Poisson statistics). 

The investigation of shot noise in a fractal conductor is 
particul arly timely in view of two different experimental 
results HH that have been reported recently. Both 
experiments measure the shot noise power in a graphene 
flake and find F < 1. A calculation [13] of the effect of 
the Pauli principle on the shot noise of undoped graphene 
predicted F = 1/3 in the absence of disorder, with a 
rapid suppression upon either p-type or n-type doping. 
This prediction is consistent with the experiment of Dan- 
neau et al. [2l|, but the experiment of DiCarlo et al. 
[2(i| gives instead an approximately do ping -independent 
F near 1/3. Computer simulations [H,[2j suggest that 
disorder in the samples of DiCarlo et al. might cause the 
difference. 

Motivated by this specific example, we study here 
the fundamental problem of shot noise due to anoma- 
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lous diffusion in a fractal conductor. While equilibrium 
thermal noise in a fractal has been studied previously 
[27L I2H , [III, it remains unknown how anomalous diffu- 
sion mi ght affect the nonequilibrium shot noise. Existing 
studies |30l . [3ll |32| of shot noise in a percolating network 
were in the regime where inelastic scattering dominates, 
leading to hopping conduction, while for diffusive con- 
duction we need predominantly elastic scattering. 

We demonstrate that anomalous diffusion affects P 
and I in such a way that the Fano factor (their ratio) be- 
comes scale independent as well as independent of df and 
a. Anomalous diffusion, therefore, produces the same 
Fano factor F — 1/3 as is known [25|, [26| for normal 
diffusion. This is a remarkable property of diffusive con- 
duction, given that hopping conduction in a percolating 
network does not produce a scale-independent Fano fac- 
tor [3(1 HH HH . Our general findings are consistent with 
the doping independence of the Fano factor in disordered 
graphene observed by DiCarlo et al. [20(. 

To arrive at these conclusions we work in the experi- 
mentally relevant regime where the temperature T is suf- 
ficiently high that the phase coherence length is <C L, and 
sufficiently low that the inelastic length is 3> L. Quan- 
tum interference effects can then be neglected, as well as 
inelastic scattering events. The Pauli principle remains 
operative if the thermal energy kT remains well below 
the Fermi energy, so that the electron gas remains de- 
generate. 

We first briefly consider the case that the anomalous 
diffusion on long length scales is preceded by normal dif- 
fusion on short length scales. This would apply, for exam- 
ple, to a percolating cluster of electron and hole puddles 
with a mean free path / which is short compared to the 
typical size a of a puddle. We can then rely on the fact 
that F — 1/3 for a conductor of any shape, provided 
that the normal diffusion equation holds locally [H, HH , 
to conclude that the transition to anomalous diffusion 
on long length scales must preserve the one-third Fano 
factor. 

This simple argument cannot be applied to the more 
typical class of fractal conductors in which the normal 
diffusion equation does not hold on short length scales. 
As representative for this class, we consider fractal lat- 
tices of sites connected by tunnel barriers. The local 
tunneling dynamics then crosses over into global anoma- 
lous diffusion, without an intermediate regime of normal 
diffusion. 

A classic example is the Sierpihski lattice [35j shown in 
Fig. [U (inset). Each site is connected to four neighbors by 
bonds that represent the tunnel barriers, with equal tun- 
nel rate T through each barrier. The fractal dimension is 
df = log 2 3 and the anomalous diffusion exponent is [6[ 
a = log 2 (5/4). The Pauli exclusion principle can be in- 
corporated as in Ref. [36[ , by demanding that each site is 
either empty or occupied by a single electron. Tunneling 
is therefore only allowed between an occupied site and 
an adjacent empty site. A current is passed through the 
lattice by connecting the lower left corner to a source (in- 
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FIG. 1: Lower panel: Electrical conduction through a 
Sierpihski lattice. This is a deterministic fractal, constructed 
by recursively removing a central triangular region from an 
equilateral triangle. The recursion level r quantifies the size 
L — 2 r a of the fractal in units of the elementary bond length 
a (the inset shows the fourth recursion). The conductance 
G — I /V (open dots, normalized by the tunneling conduc- 
tance Go of a single bond) and shot noise power P (filled 
dots, normalized by Po = 2eVGo) are calculated for a voltage 
difference V between the lower-left and lower-right corners 
of the lattice. Both quantities scale as L d i~ 2 ~ a = L l ° E2{ ' i/b) 
(solid lines on the double-logarithmic plot). The Fano fac- 
tor F = P/2eI = (P/P )(G /G) rapidly approaches 1/3, as 
shown in the upper panel. 



jecting electrons so that the site remains occupied) and 
the lower right corner to a drain (extracting electrons so 
that the site remains empty). The resulting stochastic se- 
quence of current pulses is the "tunnel exclusion process" 
of Ref. [13 • 

The statistics of the current pulses can be obtained 
exactly (albeit not in closed form) by solving a mas- 
ter equation [38[. We have calculated the first two cu- 
mulants by extending to a two-dimensional lattice the 
one-dimensional calculation of Ref. [3?J • To manage the 
added complexity of an extra dimension we found it con- 
venient to use the Hamiltonian formulation of Ref. [3t| . 
The hierarchy of linear equations that we need to solve 
in order to obtain I and P is derived in the appendix. 

The results in Fig.[T]demonstrate, firstly, that the shot 
noise power P scales as a function of the size L of the 
lattice with the same exponent df — 2 — a = log 2 (3/5) 
as the conductance; and, secondly, that the Fano factor 
F approaches 1/3 for large L. More precisely, see Fig. [2l 
we find that F — 1/3 oc L~ 15 scales to zero as a power 
law, with F — 1/3 < 10~ 4 for our largest L. 

Turning now to the application to graphene men- 
tioned in the introduction, we have repeated the cal- 
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culation of shot noise and Fano factor for the random- 
resistor network of electron and hole puddles introduced 
by Cheianov et al. jl3j |. The results, shown in Fig. [3l 
demonstrate that the shot noise power P scales with the 
same exponent L -0 97 as the conductance G (solid lines in 
the lower panel), and that the Fano factor F approaches 
1/3 for large networks (upper panel). This is a random, 
rather than a deterministic fractal, so there remains some 
statistical scatter in the data, but the deviation of F from 
1/3 for the largest lattices is still < 10~ 3 (see the circular 
data points in Fig. 



In conclusion, we have found that the universality of 
the one-third Fano factor, previously established for nor- 
mal diffusion [25l l26l . l33l . |34| , extends to anomalous dif- 
fusion as well. This universality might have been ex- 
pected with respect to the fractal dimension df (since 
the Fano factor is dimension independent), but we had 
not expected universality with respect to the anomalous 
diffusion exponent a. The experimental implication of 
the universality is that the Fano factor remains fixed at 
1 /3 as one crosses the percolation threshold in a random- 
resistor network — thereby crossing over from anomalous 
diffusion to normal diffusion. This is consistent with the 
doping-independent Fano factor measured in a graphene 
flake by DiCarlo et al. [IJ. 

A discussion with L. S. Levitov motivated this re- 
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ence Foundation NWO/FOM. We also acknowledge sup- 
port by the European Community's Marie Curie Re- 
search Training Network under contract MRTN-CT- 
2003-504574, Fundamentals of Nanoelectronics. 



FIG. 3: Same as Fig. [TJ but now for the random-resistor 
network of disordered graphene introduced by Cheianov et 
al. [l3j|. The inset shows one realization of the network for 
L/a — 10 (the data points are averaged over ~ 10 3 such 
realizations). The alternating solid and dashed lattice sites 
represent, respectively, the electron (n) and hole (p) puddles. 
Horizontal bonds (not drawn) are p-n junctions, with a neg- 
ligibly small conductance G pn ~ 0. Diagonal bonds (solid 
and dashed lines) each have the same tunnel conductance Go • 
Current flows from the left edge of the square network to the 
right edge, while the upper and lower edges are connected 
by periodic boundary conditions. This plot is for undoped 
graphene, corresponding to an equal fraction of solid (n-n) 
and dashed (p-p) bonds. 



APPENDIX A: CALCULATION OF THE FANO 
FACTOR FOR THE TUNNEL EXCLUSION 
PROCESS ON A TWO-DIMENSIONAL 
NETWORK 

Here we present the method we used to calculate 
the Fano factor for the tunnel exclusion process in the 
Sierpihski lattice and in the random-resistor network. We 
follow the master equation approach of Refs. [13, l3S| - 
The two-dimensionality of our networks requires a more 
elaborate bookkeeping, which we manage by means of 
the Hamiltonian formalism of Ref. (3f| . 



1. Counting statistics 

We consider a network of N sites, each of which is 
either empty or singly occupied. Two sites are called 
adjacent if they are directly connected by at least one 
bond. A subset S of the N sites is connected to the 
source and a subset V is connected to the drain. Each of 
the 2 N possible states of the network is reached with a 
certain probability at time t. We store these probabilities 
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in the 2 Ar -dimensional vector \P(t)). Its time evolution 
in the tunnel exclusion process is given by the master 
equation 



dt 



\P(t))=M\P(t)) 



(Al) 



where the matrix M contains the tunnel rates. The nor- 
malization condition can be written as (£|P) = 1, in 
terms of a vector (E| that has all 2^ components equal 
to 1. This vector is a left eigenstatc of M with zero 
eigenvalue 



(SI M = 0, 



(A2) 



because every column of M must sum to zero in order 
to conserve probability. The right eigenstate with zero 
eigenvalue is the stationary distribution |Poo). All other 
eigenvalues of M have a real part < 0. 

We store in the vector \P(t, Q)) the conditional proba- 
bilities that a state is reached at time t after precisely Q 
charges have entered the network from the source. Be- 
cause the source remains occupied, a charge which has 
entered the network cannot return back to the source 
but must eventually leave through the drain. One can 
therefore use Q to represent the number of transfered 
charges. The time evolution of \P(t,Q)) reads 

j t \P(t, Q)> = M \P(t, Q)> + Mi \P(t, Q - 1)) , (A3) 

where M = Mo + M\ has been decomposed into a matrix 
Mo containing all transitions by which Q does not change 
and a matrix Mi containing all transitions that increase 
Q by 1. 

The probability (£|P(i, Q)) that Q charges have been 
transferred through the network at time t represents the 
counting statistics. It describes the entire statistics of 
current fluctuations. The cumulants 



On 



d n S(t, x ) 



(A4) 



are obtained from the cumulant generating function 



S(t, X )=ln 



Q 



(A5) 



The average current and Fano factor are given by 

1= lim Ci/t, F= lim C 2 /Ci. (A6) 

t — >QC t — >00 

The cumulant generating function (jA5|) can be ex- 
pressed in terms of a Laplace transformed probability 
vector |P(t, x)> = Eo \P& Q) e* g as 



S(t,x)=ln(S|P(t,x)> 
Transformation of Eq. (|A3j) gives 
d 



dt 



\P(t,x))=M( X )\P{t,x)), 



(A7) 



(A8) 



where we have introduced the counting matrix 
Mix) =M + e x Mi. 



(A9) 



The cumulant generating function follows from 

S(t,x)=ln(S|e tM «|P(0,x)>. (A10) 

The long-time limit of interest for the Fano factor can 
be implemented as follows [38j . Let (j,(x) be the eigen- 
value of M(x) with the largest real part, and let |Poo(x)} 
be the corresponding (normalized) right eigenstate, 

M(x)\Poc(x)) = n(x)\Poa(x)), (AH) 
(S|P 00 ( X )> = 1. (A12) 

Since the largest eigenvalue of M(0) is zero, we have 

M(0)|Poo(0)) =0^/i(0) = 0. (A13) 

(Note that |Px>(0)) is the stationary distribution |P>o) 
introduced earlier.) In the limit t — > oo only the largest 
eigenvalue contributes to the cumulant generating func- 
tion, 

lim -S{t, X ) = lim ~ In [e*"W <E|Poo(x)>] = Kx)- 

t^-OQ % t—>OG t 

(A14) 

2. Construction of the counting matrix 

The construction of the counting matrix M(x) is sim- 
plified by expressing it in terms of raising and lowering 
operators, so that it resembles a Hamiltonian of quan- 
tum mechanical spins [39j . First, consider a single site 
with the basis states |0) = (J) (vacant) and |1) = (°) 
(occupied). We define, respectively, raising and lowering 
operators 





1 



1 




(A15) 



We also define the electron number operator n = s + s~ 
and the hole number operator v = 1 — n (with 1 the 2x2 
unit matrix). Each site i has such operators, denoted by 
sf , s^~, rii, and Vi. The matrix M(x) can be written in 
terms of these operators as 

M (x) = ( S K - 

(i,3) 

+ ^( e x s +-^) + ^( s r-n i ) ! (A16) 



iev 



where all tunnel rates have been set equal to unity. The 
first sum runs over all ordered pairs of adjacent 

sites. These are Hermitian contributions to the counting 
matrix. The second sum runs over sites in S connected 
to the source, and the third sum runs over sites in T> 
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connected to the drain. These are non-Hermitian contri- 
butions. 

It is easy to convince oneself that M(0) is indeed M 
of Eq. (|A1|) . since every possible tunneling event corre- 
sponds to two terms in Eq. (|A16|) : one positive non- 
diagonal term responsible for probability gain for the 
new state and one negative diagonal term responsible for 
probability loss for the old state. In accordance with Eq. 
(TM)) . the full M(x) differs from M by a factor e x at the 
terms associated with charges entering the network. 



a. Average current 

To obtain Ti we set up a system of linear equations 
starting from 

(x( X )Ti = (E| ni M( x ) \Poo(x)) ■ (A24) 

Commuting rii to the right, using the commutation rela- 
tions [lit, sf] — sf and [n», s^} = ~ s~, we find 



Extraction of the cumulants 



In view of Eq. (|Al4jl , the entire counting statistics in 
the long-time limit is determined by the largest eigen- 
value fJ,(x) of the operator (j A16|) . However, direct cal- 
culation of that eigenvalue is feasible only for very small 
networks. Our approach, following Ref. [37j |. is to derive 
the first two cumulants by solving a hierarchy of linear 
equations. 

We define 

Ti = (E| n % IPUX)) = 1 - (n n \Poo(x)) , (A17) 
U tj = U Si = (E I n inj IP,,,, ( X )> for i ? j, (A18) 
U u = 22} - 1. (A19) 

The value Ti| x= o is the average stationary occupancy of 
site i. Similarly, Uij\ x —o for i ^ j is the two-point corre- 
lator. 

We will now express /x(x) in terms of 2}. We start from 
the definition ()A1 1|) . If we act with (E| on the left-hand- 
side of Eq. (|A11[) we obtain 



(e*-l)XX 








(e x -l)E( S l 


4l^oc(x)> 






(e x ~l)E( S l 


fil-Poo(x)) 


iGS 




( e x_ 1)^(1- 


-TO- 



(A20) 



i£5 



In the second equality we have used Eq. (|A2[) [which holds 
since M = M(0)]. Acting with (E| on the the right- 
hand-side of Eq. (jAllj) we obtain just fJ,(x)> m view of 
Eq. (jA"T2|) . Hence we arrive at 



Mx) = (e x -1)E( 1 -T) 



(A21) 



ies 



From Eq. (|A21[) we obtain the average current and 
Fano factor in terms of 2} and the first derivative T- = 
dTi/dx at x = 0, 



/ = lim C x lt = //(O) = V 

t— ><X> — ' 



(i - Ti| x= o), 
Sies( 1 ~ Ti| x= o) 



(A22) 
(A23) 



Mx)T = ~ fc *T + + (e x - l) E(T - 

(A25) 

The notation X)j(i) means that the sum runs over all sites 
j adjacent to i. The number fej is the total number of 
bonds connected to site i; fc^s of these bonds connect 
site i to the source. 

In order to compute 2}| x= o we set x = in Eq. (|A25|) . 
use Eq. (|A13|) to set the left-hand-side to zero, and solve 
the resulting symmetric sparse linear system of equa- 
tions, 



ki.s — ^ ^ Tj k^Ti. 

3(i) 



(A26) 



This is the first level of the hierarchy. Substitution of the 
solution into Eq. (|A22j) gives the average current I. 



b. Fano factor 



To calculate the Fano factor via Eq. (|A23[) we also 
need T/| x=0 . We take Eq. ([£25]) , substitute Eq. ([MTj) 
for /x(x), differentiate and set x — to arrive at 

YJPk - TtT z ) - k i<s = E Tj - fciT/. (A27) 

To find C/y | x= o w e note that 

^(X)U ij = mn i n j M(x)\Poo(x)), (A28) 

and commute iij to the right. Setting x — provides the 
second level of the hierarchy of linear equations, 

= E u * + E ^ _ + fc J _ 2d «) c/ y 

+ % 5 T + fc i , 5 T :? -, i j. (A29) 

The number d,j is the number of bonds connecting sites 
i and j if they are adjacent, while dij = if they are not 
adjacent. 
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